Mon. Not. R. Astron. Soc. 000, (2002) Printed 2 February 2008 (MN I*TeX style file v2.2) 



Optical Spectroscopy of GX 339-4 - Paper I: Orbital 
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ABSTRACT 

Optical spectroscopic observations of GX 339-4 were carried out between 1998 May 
(X-ray high state) and 1999 May (X-ray low state) over 7 epochs. The equivalent 
width of Ha increased during the high state then decreased during the low state. 
The equivalent width of Hell A4686 decreased over both states. The full-width half- 
maximum of Ha, H/3, HeH A4686 and the Bowen blend increased from 1998 to 1999 
indicating that the emission line regions moved closer to the compact object. H/? shows 
a redshifted absorption feature at A4880 at all epochs. This line remains unidentified. 
Analysis of individual spectra from 1998 May 28-31 show modulation of the radial 
velocities, equivalent width and V/R ratios of Ha on the 14.86 hour orbital period. 
The equivalent width of HeH A4686 also varies on the orbital period. This is the first 
time since the study of Callanan ct al. ( 1992) that spectroscopic data has confirmed 
the orbital period. The semi-amplitude of the Ha radial velocities is Ki = 14 km s^^. 
Hence the mass function — 2 x 10^"^ Mr. 
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1 INTRODUCTION 

GX 339-4 was di scovered by 05*0 - 7 in 1973 
iMarkert at all Il973|). The secondary star has a mag- 
nitude B > 21 Jkutchings Cowlev fc CramptonI Il98lt 
lllovaiskv fc ChevalieJ ll98J) . Consequently the secondary 
has never been directly observed and its spectral type never 
identified. 

An orbital period of 14.86 hours (0 .61916±0.0027 days) 
was derived by ICaUanan et"ai] (Il99l hereafter CCHT92) 
from R- and 7-band photometry taken during an X-ray ofi' 
state. Two sets of spectra obt a ined in the high state by 
ICowlev Crampton fc Hutchingj (I^S^) were folded on this 
period by CCHT92. The first data set did not show any 
modulation but the second set did, giving a semi-amplitude 
Ki = 78 ± 13 km s~^ and systemic velocity 7 = -62 ± km 
s~^. The mass of the compact object was estimated to be 
Ml ^ 1 Mq based on these measurements and assuming 
the binary inclination i ^ 70°. This estimate of Mi is well 
below t he theoretical maximum mass of a neutron st ar ('^ 
3.2 M^. IChitre fc H artic 1976: F riedman fc Ipserll987l) . GX 
339-4 has been classified as a black hole X-ray binary based 
purely on X-ray characteristics and behaviour which are very 



* Now at Astronomy Department, Yale University. E-mail: bux- 
ton@astro.yale.edu 



similar to another black- hole candidate Cyg X-1. If GX 339- 
4 is not black hole then this places the method of black-hole 
classification using X-ray data into serious question. 

CCHT92 proposed that the optical modulation is due 
to eclipses of an X-ray heated companion and/or accretion 
disk as it was seen in both the off and high state data. This 
would require i ^ 60° which has not been supported by 
other optical or X-ray observations. 

The orbital period has not been confirmed by either op- 
tical spectroscopy or photometry since the CCHT92 study. 
This means that the orbital period is wrong, that the orbital 
period is only detectable under certain physical conditions 
in the binary or that the data quality of subsequent studies 
has not been suflicient. 

The main aim of this study was to determine the orbital 
period of GX 339-4 from emission lines in optical spectra. 
This method is notoriously diflicult. As the emission lines 
originate from the accretion disk information such as the 
radial velocities are more difiicult to extract since their am- 
plitudes are much smaller and they may be contaminated 
by other fiux sources. Since we have no spectral signature of 
the secondary star we are left with no other choice but to 
use disk emission lines to derive the orbital period. 

In this paper we present the results of analysis con- 
ducted on the emission lines in GX 339-4. We discuss the 
spectral changes seen over the epochs and how these changes 
relate to the X-ray behaviour over the same period. In ad- 
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dition we present results of analysis aimed at detecting the 
orbital period from disk emission line variations. 



2 OBSERVATIONS AND DATA REDUCTION 

A total of 51 low dispersion and 59 high dispersion spec- 
tra were obtained between 1998 May and 1999 May. An 
observation log is given in Tabled Observations were per- 
formed using the Double-Beam Spectrograph (red and blue 
arms) on the Australian National University 2.3m telescope 
and the Royal Greenwich Observatory spectrograph on the 
Anglo- Australian Telescope. 

Figure Q shows the epochs of observations with respect 
to the soft (1.5-12 keV) and hard (20-100 keV) X-ray obser- 
vations. Our observations of GX 339-4 were obtained during 
high and low X-ray states. 

Data reduction was performed using TRAP. Each spec- 
trum was corrected for electronic bias and pixel-to-pixel vari- 
ations. The background was fitted using a low-order polyno- 
mial. Wavelength calibration was performed via FeAr (blue) 
and NeAr (red) arc spectra from which line positions were 
fitted with a fourth-order polynomial. The largest rms error 
on the wavelength calibration for blue spectra was 0.01 A 
(^ 1 km s~^) in high resolution spectra and 0.8 A (^ 50 km 
s~^) in low resolution spectra. In red spectra the largest rms 
error was 0.01 A (~ 0.5 km s~^) in high resolution spectra 
and 0.4 A (~ 18 km s~^) in low resolution spectra. Telluric 
lines were removed using a smooth-spectrum standard star 
applying the method outlined in BcsscU (1999). For nights 
which had clear observing conditions spectra were fiux cal- 
ibrated using the standard stars Feige 110, LTT 7379, HD 
49798 and HD 127493. 



3 SPECTRA IN HIGH AND LOW X-RAY 
STATES 

Figure |5| shows combined, normalised blue and red spectra 
from the high and low X-ray states. Only low dispersion 
spectra are shown here to present full wavelength coverage. 
Spectra were combined with weighting on individual spectra 
according to their signal-to- noise ratio (SNR). 

Many emission lines are present in the high state. The 
strongest lines are Hq and Hell A4686. Other emission lines 
visible include the H7, H/3, the Bowen blend (~ A4640 A), 
Hell A5411, Hel AA4922, 5875 (blended with NaD), 6678 
and 7281 and the Paschen series. 

In the low state most of the lines are still visible. The 
most dramatic flux decrease from the high state to the low 
state is seen in H5, the Bowen blend, Hell A4686 and the 
lower members of the Paschen series. 

The absorption lines present are Call-K (A3934), the 
NaD AA5889, 5895 lines and the diffuse lines AA4430 and 
6283. These are of interstellar origin. 

FigureOshows an expanded view of the dominant emis- 
sion lines in combined spectra over the high and low states. 
The combined spectra in this case include both high and 
low resolution data with weighting on the SNR of individ- 
ual spectra. 

An increase in line width is most obvious in Ha but less 
so in other lines. The red wing of Hell A4686 and H/3 have 
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Figure 3. Combined high and low resolution spectra from high 
(histogram) and low states (line) expanded around the most dom- 
inant emission lines. 



slightly increased. The Balmer and Hel lines are asymmet- 
ric on the blue side in the high state. They become more 
symmetrical in the low state. 

The H/3 emission line shows a redshifted absorption dip 
in both the high and low state data. Further discussion on 
this feature is given in H4.5I 

The similarities in flux behaviour between the Bowen 
blend and Hell A4686 suggest that both emission lines 
arise from the Bowen Flourescenc e Mechanism (BFM, 
iMcChntock Canizares fc TarteJll975ll . Other emission lines 
expected from the BFM are OIII AA3133 and 3444. Spectra 
are too noisy in this region to ascertain whether these lines 
are present. Weaker emission from NIII A4097-4103 is also 
expected but would be blended with H(5. This could explain 
why H5 shows the most flux decrease of the Balmer lines. 
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Table 1. Observation log of GX 339-4 optical spectroscopy. 



Date 


Telescope 


Grating 


Resolution 


Exposure 


Wavelength Range 


No. of 






(1/mm) 


(A) 


(sec) 


(A) 


Spectra 


1998 May 26,27 


ANU 2.3m 


300/316 


4.4 


900,1800,2700,3600 


3300-9500 


20 


1998 May 28-31 


ANU 2.3m 


1200 


2.2 


3600 


4600-5300, 6100-7050 


25 


1998 Aug 20-23 


ANU 2.3m 


1200 


2.2 


2000 


4400-5100, 6200-7150 


17 


1999 Mar 11 


ANU 2.3m 


600 


4.4 


900 


3150-8800 


5 


1999 Mar 24-27 


ANU 2.3m 


300/316 


4.4-8.4 


900 


3300-9450 


20 


1999 Apr 13 


AAT 3.9m 


1200 


1.5 


1800 


6480-7290 


3 


1999 Apr 17, 18 


ANU 2.3m 


1200 


2.2 


1800 


4250-4970, 6050-7000 


14 


1999 May 11 


ANU 2.3m 


300/316 


4.4 


1800 


3300-9450 


6 




Figure 1. Top: The ASM day-averaged count rate (counts/sec) of GX 339-4 in 1.5 - 12 keV band (1 Crab ~ 75 counts/sec). Quick-look 
results provided by the ASM/RXTE team. Bottom: BATSE photon fluxes in 20 - 100 keV band. Arrows indicate epochs of optical 
observations. 'Low', 'High' an d 'Off' refer to the X -ray states. Transition times (indicated by the vertical dashed lines) and X-ray state 
classifications were taken from lBelloni et all 1 19991) and lKong et al.l l2000l) . 



The decrease in Hell A4686 and Bowen blend emission 
is most likely due to fewer soft X-rays available for pho- 
toionisation. There must be enough, however, to sufficiently 
excite the lower ionisation lines of Hel and the Balmer series. 
Parts of the disk which were heated in the high state may be 



allowed to cool sufficiently in the low state for Balmer and 
Hel emission to occur closer to the compact object. This 
would explain the increase in line widths (see also i|4.41 . 

In the next sections we will present quantitatively how 
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Figure 2. Left three panels are blue spectra and right three panels red spectra of GX 339-4 in the high (solid line) and low (histogram) 
X-ray states. Visible lines are indicated. Low resolution data only are shown here to present full wavelength coverage. 'BB' is the Bowen 
blend. 'Telluric' indicates region of strongest telluric lines which may not have been removed perfectly. 



the emission lines vary on long (days to months) and short 
(hours) timescales. 

The analysis of the spectra has been divided into two 
parts. The first deals with combined spectra from which the 
distance, extinction, magnitudes and colours are measured. 
Consideration is also given as to how the equivalent width 
(EW) and full-width half maximum (FWHM) of the emis- 
sion lines vary in comparison to the soft X-ray behaviour 
between 1998 and 1999. 

The second part of the analysis is concerned with indi- 
vidual spectra. The analysis aims to determine any orbital 
modulation which may be present in the properties of the 
emission lines. 

In a companion paper (Paper II) we will explore further 
the line profile behaviour and how the variations may be 
related to various sources in the binary. 



4 ANALYSIS OF COMBINED SPECTRA 
4.1 Distance and Extinction 

Previous measurements of the distance D and extinction 
E{B - V) to GX 339-4 are summarised in Tabled 

Spectra were combined for nights when the NaD and 
Call-K lines were included in the spectral range. Four or five 
measurements of the EWs of NaD and Call-K were made 
from the combined spectra for each night. Measurements 
were obtained by fitting single-peaked Gaussians to the line 
profiles using splot in IRAF, varying the placement of the 
continuum. They are listed in Table |3 Errors were derived 
from the scatter in the measurements. 

The grand average EW of the NaD fines is 4.2 ± 1.0 
A and for Call-K 1.6 ± 0.3 A. The former result agrees 



with Soria et alj (|19 99[) bu t is higher than ICorbet et alJ 
Jl987) and iGrindlavl lll97Sl) . The latter agrees weU with 



ICowlev et alJ ^ll987^ . 

Using the method described in lAllerJ lll97 3l) we find that 
73 ~ 4 ± 1 kpc. For Ca l l-K w e used the empirical relation- 
ship from lBeals fc Okd lll953l) giving D ~ 4 ± 1 kpc. This 
agrees with the result from the NaD lines. Beals & Oke also 
provide a relationship for the NaD lines. In this case D ~ 
7 ± 2 kpc which is larger than the other results but just 
consistent within the errors. Since the first two results agree 
well with one another and with previous studies we have 
adopted D = 4 ± 1 kpc. The other measurement should 
not be ignored, however, and indeed highlights the inconsis- 
tency which various methods of distance determination via 
interstellar lines may give. 

E( B — V) was calculat ed using the empirical relation- 
ship of iBarbon et alJ ([1990') and was found to be ~ 1.1 ± 
0.2. This is consistent with Grin dlay (19?.d)- Using the in - 
terstellar extinction relationship of lSavage fc Mathij il979l) , 
Av = 3.1E{B -V)= 3.3 ± 0.7. 



4.2 Continuum Variations 

Combined, fiux-calibrated spectra are shown in Figure 0| 
Individual spectra were averaged with weighting according 
to their SNR. The spectra were dereddened using Ay = 

3.3 ±0.7 (S3J- 

The continuum is blue which is typical of XRBs during 
outburst. The continuum flux decreased steadily over the 
year of observations. The continuum reddened as the binary 
went from the high state to low state (see also H4.3II . This is 
due to less optical flux originating from the accretion disk as 
fewer X-ray photons are available to be reprocessed by the 



Optical Spectroscopy of GX 339-4 - Paper I 5 



Table 2. Previous measurements of distance and extinction of GX 339-4. 



Reference 


Distance 


EW oi NaD 


EW ol Call-K 


E{B — V) 




(kpc) 


(A) 


(A) 


Doxev et al. (1979) 


~ 4 








Grindlav (1979) 


~ 8 


2.2 ± 0.4 




~ 1.2 


Mauche & Gorenstein (1986) 


1.33 








Makishima et al. ("19861 


3.5 








Ilovaiskv et al. ( 19861 








0.7 ± 0.1 


Cowlev et al. ("19871 


3.5-4 




1.3 




Corbet et al. (19871 




1.0 






Predehl et al. (19911 


1.3 








Zdziarski et al. (19981 


4 ± 1 








Soria Wu & Johnston (19991 




~ 3.7 






Shahbaz Fender & Charles (20011 


5.6 









Table 3. Measured equivalent widths of NaD and Call-K inter- 
stellar absorption lines from GX 339-4 spectra. 



Date EW of NaD EW of Call-K 

(A) (A) 



1998 May 26 


3.8 


± 


0.4 


1.6 


± 


0.4 


May 27 


5.5 




1.5 


1.7 




0.4 


1999 Mar 11 


4.0 




0.6 








Mar 24 


3.8 


± 


1.2 








Mar 25 


4.3 


± 


1.1 


1.5 




0.3 


Mar 26 


4.5 


± 


0.8 


1.6 


± 


0.2 


May 11 


3.7 




0.5 








Average 


4.2 


± 


1.0 


1.6 


± 


0.3 
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Figure 4. Flux-calibrated spectra over some observed epochs. 
Flux units are ergs cm~^ s~^ A~^. Spectra have been heavily 
smoothed to accentuate the continuum behaviour. The continuum 
decreases and reddens from 1998 to 1999 but the slope can vary 
greatly from night to night as seen with 1999 March 25 and 26 
spectra. 



Table 5. Colours of GX 339-4 measured from combined, dered- 
dened, flux-calibrated spectra. 



Date 






V) 


iv 




R) 


(V 




I) 


1998 May 26 


-0.23 


± 


0.03 


0.27 




0.05 


0.10 


± 


0.10 


1999 Mar 11 








0.44 


± 


0.05 








Mar 25 


0.02 


± 


0.03 


0.41 


± 


0.05 


0.44 


± 


0.10 


Mar 26 


-0.15 


± 


0.03 


0.35 


± 


0.05 


0.31 


± 


0.10 


May 11 


0.06 


± 


0.03 


0.45 


± 


0.05 


0.56 


± 


0.10 



outer disk. From night to night, however, the slope varied 
dramatically as seen with spectra from 1999 March 25 and 
26. On these nights the slit width was 2 arcseconds and 
seeing conditions were noted as being 1.5 - 2 arcseconds. 
Hence most of the flux from the system would have been 
recorded and changes in the continuum are probably real. 



4.3 Optical Magnitudes and Colours 

B, V, R and I magnitudes and {B-V), (V-R) and (V-I) 
colours were derived from the combined fiux-calibrated, 
dereddened spectra using sbands in IRAF. The results are 
shown in Tables |1| and |^ Errors in the magnitudes were 
determined by finding the scatter between the sensitivity 
functions of the standard stars used in the flux calibration. 
The errors in the colours were derived by measuring the 
difference in slopes between the sensitivity functions. 

The trend from 1998 to 1999 suggests that the magni- 
tudes have decreased. The errors, however, are larger than 
the greatest difference between the magnitudes. Hence we 
cannot say whether this trend is real. 

All three colour indices increased from 1998 May to 1999 
May (see Figure |SJ. In this case the errors are smaller than 
the greatest difference between 1998 and 1999 suggesting 
that this trend is real. There is, however, large changes be- 
tween 1999 March 25 and 26. Therefore we cannot deter- 
mine exactly how much the colours change between 1998 
and 1999. 
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Table 4. Magnitudes of GX 339-4 measured from combined, dereddened, flux-calibrated spectra. 



Date 




B 






V 






R 










1998 May 26 


12.89 


± 


0.75 


13.12 


± 


0.75 


12.85 


± 


0.75 


13.02 


± 


0.75 


1999 Mar 11 








13.35 


± 


0.75 


12.91 


± 


0.75 








Mar 25 


13.51 


± 


0.75 


13.49 




0.75 


13.08 


± 


0.75 


13.05 




0.75 


Mar 26 


13.57 


± 


0.75 


13.72 




0.75 


13.37 


± 


0.75 


13.41 


± 


0.75 


May 11 


14.13 


± 


0.75 


14.07 


± 


0.75 


13.62 


± 


0.75 


13.51 


± 


0.75 
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Figure 5. Optical colours measured from combined flux- 
calibrated, dereddened spectra. Solid boxes are {B — V), open 
boxes {V — R) and open triangles (V — I). 

4.4 Emission Line Variations over Epochs 

Spectra from each epoch were combined and normalised to 
the continuum. Low and high resolution data were separated 
and remain so for the rest of this analysis. 

The EW and FWHM of the strongest emission lines 
were measured by fitting single-Gaussian profiles using splot 
in IRAF. For each spectrum three or four measurements 
were made of each line varying the placement of the con- 
tinuum. The average EW and FWHM were calculated from 
these measurements and the errors determined from their 
scatter. 



4-4- i Equivalent Width 

The EW of Hq, H/3, Hell A4686 and the Bowen blend for 
each epoch are listed in Table |S| and are plotted in Figure |S| 
together with the soft X-ray flux. 

The 1998 May data show a decrease of EW in Ha, H/? 
and HeH A4686 from low resolution to high resolution data 
by ~ 0.7, 0.2 and 0.6 A, respectively. There is an increase 
in the Bowen blend EW but the errors are quite large. 

The EW of Ha increased in 1998 August by - 1.8 A 
while that of Hell A4686 and the Bowen blend decreased by 
0.8 and 0.6 A, respectively. 



Figure 6. Equivalent widths (A) of the strongest emission lines 
in GX 339-4 from combined spectra of observed epochs. Bottom 
panel shows soft X-rays for comparison. 'High' and 'Low' refer to 
the X-ray states. 

In the low state one clearly sees a decreasing trend of 
EW in Ha, H/3 and HeH A4686 and an increase in the Bowen 
blend, although the Bowen blend again has large errors. 

We have seen that changes in the continuum slope are 
possible (see Figure 2J- K the slope of the continuum be- 
comes more negative (that is, bluer) then the continuum 
flux may decrease the equivalent width of HeH A4686 and 
increase that of Ha. Hence the continuum may vary suffi- 
ciently to account for some of the changes seen here. 

The decrease in EW during the low state suggests that 
the line fiux must be decreasing faster than the continuum. 

Without simultaneous photometry it cannot be said 
which of the continuum or line fiux variations are the major 
causes to the EW changes over the year. Analysis of future 
spectroscopic observations will benefit greatly if obtained si- 
multaneously with photometry so that we may discern the 
cause of the EW changes. 

4.4.2 Full- Width Half Maximum 

The FWHM of the strongest emission lines are given in Table 
Qand plotted in Figure |7] with the soft X-ray fiux. 

The FWHM decreases between the low and high reso- 
lution data of 1998 May in Ha, H/3 and Hell A4686 by ~ 



Optical Spectroscopy of GX 339-4 - Paper I 7 

Table 6. Equivalent widths (A) of the strongest emission lines in GX 339-4 from combined spectra of observed epochs. 



Epoch 


HJD 


Ha 


H/3 


Hell A4686 


Bowcn Blend 


1998 May (low res.) 


2450960.61940 


6.8 ± 0.2 


2.8 ± 0.2 


6.3 ± 0.2 


3.3 ± 0.2 


1998 May (high res.) 


2450963.08762 


6.2 ± 0.2 


2.6 ± 0.1 


5.7 ± 0.1 


3.5 ± 0.2 


1998 August 


2451047.47701 


8.0 ± 0.3 


3.0 ± 0.8 


4.9 ± 0.2 


2.9 ± 0.3 


1999 March 


2451257.19183 


8.5 ± 0.1 


3.1 ± 0.2 


5.9 ± 0.1 


3.0 ± 0.2 


1999 April 


2451284.61028 


7.6 ± 0.3 


3.3 ± 0.5 


5.1 ± 0.3 


3.3 ± 0.3 


1999 May 


2451310.24667 


7.4 ± 0.4 


2.4 ± 0.3 


3.8 ± 0.4 





Table 7. Full-width half maxima (km s ^) of the strongest emission lines in GX 339-4 from combined spectra of observed epochs. 



Epoch 


HJD 


Ha 


Uf3 


HeH A4686 


Bowen Blend 


1998 May (low res.) 


2450960.61940 


667 ± 14 


796 ± 31 


755 ± 13 


1130 ± 51 


1998 May (high res.) 


2450963.08762 


564 ± 8 


748 ± 19 


717 ± 6 


1210 ± 45 


1998 August 


2451047.47701 


507 ± 9 


592 ± 120 


672 ± 13 


976 ± 71 


1999 March 


2451257.19183 


709 ± 5 


808 ± 25 


868 ± 12 


1400 ± 52 


1999 April 


2451284.61028 


667 ± 14 


833 ± 68 


851 ± 26 


1750 ± 100 


1999 May 


2451310.24667 


699 ± 23 


826 ± 56 


883 ± 64 
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Figure 7. Full-width half maxima (A) of the strongest emission 
lines in GX 339-4 from combined spectra of observed epochs. 
Bottom panel shows soft X-rays for comparison. 



100, 140 and 45 km s~\ respectively. In all lines the FWHM 
decreased from 1998 May to 1998 August then increased in 
1999 by as much as ~ 190 km s~^ for Ha. 

The increase in FWHM from the high state to the low 
state is most noticeable in the line wings (Figure 01 . This 
suggests that the line emitting regions are moving closer to 
the compact object from the high state to low state. As 
stated above (3SJ, the disk may be cooling in the inner disk 
regions during the low state as there are fewer soft X-rays 
available to heat the disk. Therefore regions closer to the 
compact object may emit Balmer lines durin g thi s time. 

It is worth noting that lSoria et alJ lll999fl and lWu et alJ 



i200ll) saw no change in the line width of Ha between 
the high and low X-ray states but that of Hell A4686 in- 
creased with X-ray hardness. In these observations a change 
in FWHM is detected in all of the prominent emission lines, 
including Ha, from the high to low state. It would seem, in 
this case at least, that all the strongest emission lines are 
correlated with X-ray hardness. 

The Bowen blend FWHM is substantially larger than 
Hell A4686. This would imply that the Bowen blend emis- 
sion originates closer to the compact object. The BFM, how- 
ever, requires that Hell A4686 and the Bowen blend emission 
originate from the same place since some of the Hell La pho- 
tons are absorbed and remitted by OIII and NHL Since it 
is a blend of several lines, the FWHM of the Bowen blend 
may be 'artificially' increased. 

4.5 H/3 Absorption Feature 

An interesting feature in the combined blue spectra is the 
redshifted absorption dip at H/3 which is seen at all epochs. 
No other emission lines exhibit such absorption. Figure |21 
shows combined low-dispersion spectra. As can be seen. Ha 
has a redshifted absorption during the low state. When low 
and high dispersion spectra are combined (Figure El the 
absorption dip is not evident. Hence the absorption feature 
near Ha could be an artifact of combining low dispersion 
spectra only whereas the redshifted absorption component 
is present near H/3 in all spectra. 

Figure |H| shows the combined spectra from all epochs 
expanded around H/3. Several Gaussian fits were made to 
the H/3 emission and absorption lines. An example fit is 
also shown in Figure |H| The Gaussian fits have the following 
properties: H/3 Gaussian: FWHM = 780 ± 200 km s"\ EW 
= 2.67 ± 1.0 A, Ac = 4861.94 ± 0.03 A; absorption Gaus- 
sian: FWHM = 1200 ± 280 km s-\ EW = 0.90 ± 0.03 A, 
Ac = 4880.49 ± 0.14 A. The errors shown were derived from 
the scatter in the measurements. 

The same feature has been observed in V404 Cyg by 
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Figure 8. Redshifted absorption near H/3. The solid line is the 
combined spectrum from all epochs centred on H/3. The emission 
line to the right is Hel A4922. The dot-dash line is the Gaussian 
fit to H/3 emission line. The dashed line is the Gaussian fit to 
absorption trough. The combined Gaussian is shown by the long- 
dashed line. 



ICasares et al.l il993l) who attribute it to a metallic blend of 
CrI, Nil and Fel at A4889. Casares et al. found the spectral 
type of the secondary in V404 Cyg to be G9±l V to K0±4 
III derived by comparison of Cal and Pel lines in the range 
AA 6350-6530 to template star spectra. No absorption lines 
were seen in this spectral range for GX 339-4. Thus it seems 
unlikely that the observed absorption near H/3 originates 
from the secondary star. 

If the absorption line is an inverted P-Cygni profile then 
we would expect to also see it in Ha. As Ha does not con- 
vincingly show a redshifted absorption feature in all spectra 
this explanation is also doubtful. 

Future observations should endeavour to see if this fea- 
ture is still present near H/3 and whether Ha exhibits the 
same or otherwise. Obtaining higher SNR ratio data may 
enable us to ascertain the source of this absorption. 
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Figure 9. Individual high resolution blue (top) and red (bottom) 
spectra of GX 339-4 on 1998 May 28 and 29. The emission lines 
visible from left to right are, in the blue spectra, Bowen blend. 
Hell A4686 and H/3 and, in the red spectra. Ha and Hel A6678. 



5 ANALYSIS OF INDIVIDUAL SPECTRA 

An example of low and high resolution individual spectra 
are plotted in Figures 151 and [Till Inspection of the spectra by 
eye show that the line profiles vary on a timescale of hours. 
One can see in particular the peak flux in Hell A4686 vary 
on the blue side on May 29. The Hell A4686 profile on May 
28 looks more symmetrical and the double-peaks are not as 
clear. In Ha it is difficult to discern two well-discerned peaks 
at any time although one may argue that they are present 
in the first few spectra on May 29 and perhaps May 31. Ha 
is more round-topped and symmetric on May 28 then be- 
comes asymmetric on the blue side on May 29. The spectra 
repeat this profile behaviour the following two nights (Fig- 
ure Since the profile characteristics repeat themselves 



this strongly suggests that the profile variations are related 
to the orbital period. 

The following analysis attempted to extract quantita- 
tive information on the changes in the line profiles to see if 
they indeed corresponded to an orbital period and whether 
this agreed with the published orbital period of CCHT92. 

At this point we note that the profile shape changes 
from 1998 to 1999 when we begin to see red-sided asym- 
metry. This is studied in more detail in Paper II. Plots of 
indivi dual spectra from other epochs can be found iii iBuxtonl 




The analysis was conducted in three ways. 
First a single Gaussian profile was fitted to the most 
prominent emission lines using splot in IRAF. Several mea- 
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Another way of quantifying the hne profile variations is 
by taking the ratio of fluxes from the blue (or violet) and 
red sides of the emission line, also known as the V/R ratio. 
The advantage of this method is that no assumptions are 
made on the shape of the line profile. The V/R ratio was 
calculated for Ha and Hell A4686. The dividing line between 
the blue and red fluxes was taken to be the rest wavelength, 
that is 6562.80 A for Ha and 4685.68 A for HeH A4686. 

A period search was conducted on the EW, FWHM, 
radial velocities and V/R ratios. The search was performed 
over a restricted range of frequencies to eliminate possible 
aliases. The frequencies searched over ranged from three 
times the exposure of the individual spectra to one third 
of the epoch coverage. 

Due to the size of the dataset it is not possible to show 
every single result of the analysis. Hence, we will first discuss 
briefly the datasets for which the results are not presented 
here. 

In only a few cases, epochs did not have sufficient tem- 
poral coverage such that the published orbital period was 
not sampled in the periodogram. This was always true for 
the low resolution f998 May dataset. 

The SNR of individual spectra at many epochs was 
quite low. The SNR of high resolution spectra in 1998 May 
was 50-60 compared to 10-20 for high resolution spectra in 
1998 August. The errors in EW, FWHM and radial veloci- 
ties were too large to give any meaningful results in the low 
SNR cases. 

Analysis of the double-gaussian fits did not show any 
modulation in the radial velocities, EW, FWHM or V/R 
ratios. 

Most epochs which had a SNR of ~ 40 did not have a 
sufficient number of spectra to cover all orbital phases. 

Hence, the discussion of results in the following sections 
concentrate on those from the high resolution data of 1998 
May 28-31. This dataset provided the best temporal cov- 
erage, had the highest SNR and had the most number of 
spectra. 

For results on analysis of other epochs the reader is 
referred to .Buxton. ^2002) ■ 



Figure 10. As for Figure|2lfor 1998 May 30 and 31. 

surements were taken of each line in each spectrum varying 
the placement of the continuum. The average of the EW, 
FWHM and line centres (leading to radial velocities) were 
calculated from these measurements. Single-Gaussian pro- 
files do not perfectly fit the line profiles, particularly in the 
central regions of the emission line. As the profiles fitted the 
line wings well we were satisfied that the FWHM and line 
centres were measured accurately. To test how well the EW 
was measured, we used the "e" option in splot within IRAF 
and found that there were virtually no differences between 
the two methods. 

The second method of analysis used double-Gaussians 
to fit profiles in cases where the two peaks were well dis- 
cerned. Several measurements were made varying the con- 
tinuum and Gaussian centre placements. The average EW, 
FWHM and radial velocities of both Gaussian components 
were calculated. 



5.1 Radial Velocities 

Periodograms created from the Ha, HeH A4686 and the 
Bowen blend radial velocities measured via single-Gaussian 
profiles are shown in Figures [TTl to [HI The 68% and 99% 
confidence levels shown were calculated using the method 
outlined in lLampton Margon fc Bowverl il976ft for the two- 
parameter model case. 

There is a strong peak in all periodograms at the 99% 
confidence level. The peak for Hq falls at the published or- 
bital period. The periodogram of Hell A4686, however, has a 
broad double-peak around the published orbital period. The 
Bowen blend shows a strong peak at the published orbital 
period. There are, however, peaks at other periods which 
are as significant. 

Before folding the radial velocities on the orbital pe- 
riod it was first necessary to calculate a new ephemeris. The 
ephemeris given by CCHT92 is Yorb = 0.61916±0.0027 days 
and To (JD, minimum light) = 2447000.539 ± 0.015. The 
error on the period is large enough to render the ephemeris 
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Figure 11. Top: Periodogram of Ha radial velocities measured 
from single-Gaussian line profile fits. Data is high resolution spec- 
tra from 1998 May 28-31. Arrow shows where the published pe- 
riod of 14.86 hours lies, y-axis is the normalised, reduced ■ 
Long dashed line shows 99% confidence level level and dotted 
line the 68% confidence level. Bottom: Radial velocities of H« 
measured from single-Gaussian profile fits phased on the 14.86 
hours period. Points were grouped into 0.10 phase bins. Solid 
curve is best-fit sinusoid derived from the least-squares sine fit 
to the binned data. Errorbars are the standard deviation derived 
from the least-squares fit. 

useless after ~ 140 days and, therefore, could not be used 
to correctly phase this dataset. 

A new ephemeris was determined from the radial veloc- 
ities of Ha which were measured from single-peaked Gaus- 
sians. To obtain the phase shift a period search was con- 
ducted on the data with the period restricted to the pub- 
lished period of 0.61916 days. The amount of time corre- 
sponding to this phase shift was subtracted from the time 



Figure 12. As for Figure E] for Hell A4686. 

of the first observation, giving a new To (JD) = 2450961.8773 
± 0.0078. The error in To was calculated by conducting the 
period search again on the period with the error subtracted 
and added and finding To in both cases. The error given is 
the larger of the two. 

The orbital period could not be refined as the errors on 
the radial velocities over the various epochs were too large. 

Phase 0.0 corresponds to the time when the radial ve- 
locities cross from the blue to the red. Since the majority of 
the emission line fiux originates from t he accretion disk as 
indic ated by the double-peaked profiles teuandll972HSmakl 
Il969l) the radial velocity variations most likely relate to the 
disk. Therefore, phase 0.0 corresponds to the inferior con- 
junction of the primary. One must use caution, however, as 
line asymmetries caused by, say, a hotspot may shift the line 
centroid and give a false indication of the time of maximum 
and minimum of the radial velocities. This would require 
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Table 8. Parameters of best-fit sinusoid curves to radial velocities 
of emission lines. 



Line Ai 7 

(km s^^) (km s^^) 

Hq 14 ± 3 -18 ± 2 

Hell A4686 20 ± 20 8 ± 10 

Bowen blend 30 ± 20 -10 ± 20 
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Figure 13. As for Figure ITT1 for the Bowen blend. 

more detailed analysis involving line profile modelling which 
is pursued in Paper II. For now we take phase 0.0 to be the 
point at which the accretion disk is between the observer 
and the secondary star. 

The radial velocities of Ha, Hell A4686 and the Bowen 
blend were folded on the published orbital period using the 
new To. The results are also shown in Figures 1111 to 1131 
Points were binned by 0.10 in phase. The binned data were 
subjected to a least-squares sine fit. The solid curves in the 
figures show the best-fit sinusoid. Their parameters are listed 
in Table |H] The errors on the semi-amplitudes and systemic 
velocities were obtained from the fit. The errorbars on the 
radial velocities are the standard deviations obtained from 
the best-fit sinusoids. 

The Ha radial velocity curve shows clear modulation on 
the published orbital period of 14.86 hours. The maximum 
is at phase 0.25 and the minimum at phase 0.75. 



The radial velocity curve of Hell A4686 is not as well 
defined as Hq and is obviously why the periodogram did 
not show an isolated peak at the published orbital period. 
Consequently the best-fit sine curve does not describe the 
data well. 

The Bowen blend shows a much cleaner radial velocity 
curve than Hell A4686. The maximum, however, is at phase 
0.5 and minimum at phase 0.0. This is 0.25 out of phase with 
respect to the Ha radial velocities. In addition, the velocity 
amplitude is much larger than Ha and Hell A4686. Since 
the Bowen blend is a blend of several lines which may vary 
in strength with time it is not a reliable velocity tracer. 

Previous studies which use emission lines to 
determine 7 show that it is sensitive to distor- 
tions in the line profile and cannot be relied upon 
JOrosz. Bailvn. Remillard. McClintock fc Foit3 Il994l : 
iGarcia. Callanan. McClintock fc Zhaol Il99d) . Since the 
lines are obviously asymmetric the 7 derived here will not 
be discussed further but have been shown in Table |H| for 
completeness. 

CCHT92 derived Ki = 78 ± 13 km s~^ f rom t he Ha 
radial velocity data obtained bv lCowlev et al.l 1^83). This 
value is much higher than that obtained here. A lower radial 
velocity amplitude decreases the mass function and increases 
the upper limit to the compact object mass. Using Ki = 14 
km s~^ from Ha, and Porh = 14.86 hours, the mass function 
is /(M) = 2 X 10"* M©. Assuming M2 = 0.7 Mp, (K5 V 
star) and i ^ 70°, Mi ^ 40 Mq. If i = 15° JWu et alJl26ol 
then Ml = 5 M© . This would imply that the compact object 
in GX 339-4 is a black hole. One must use caution, however, 
as previous studies have shown th at using em ission lines to 
measure Ki may not be accurate iWaddll985l) . In addition, 
the binary inclination is not well constrained although it is 
likely to be low. 

5.2 Equivalent Width 

Periodograms of the EW for Ha and Hell A4686 are shown 
in Figures 1141 and 1151 Both cases show peaks at the 99% 
confidence level at the published orbital period. The Bowen 
blend periodogram did not show any significant peaks (see 
iBuxton 2002). 

The EW of Ha and Hell A4686 were folded on the or- 
bital period using To calculated in tl5.1l Points were binned 
by 0.10 in phase. Figures |14I and |l5l convincinglv show a 
modulation of the EWs on the orbital period for both lines. 
The sinusoid curves shown in the figures are the best fit to 
the data found from a least-squares sine fit. Their parame- 
ters and errors are given in Table |5] The errors on the EW 
are the standard deviations calculated from the least-squares 
sine fit. 
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Figure 14. Top: Periodogram of Ha EW measured from single- 
Gaussian profile fits. Data is high resolution spectra from 1998 
May 28-31. Arrow shows where the published period of 14.86 
hours lies. j/-axis is the normalised, reduced ■ Long dashed line 
shows 99% confidence level level and dotted line the 68% confi- 
dence level. Bottom: EW of Ho folded on the published period. 
Data points were binned by 0.10 in phase. Solid curve is the best 
fit found from the least-squares sine fit to binned data. Errorbars 
are the standard deviation derived from the least-squares fit. 



Table 9. Parameters of best-fit sinusoid curves to the EW of Ha 
and HeH A4686. 



Line 


Semi-amplitude 


Mean 




(A) 


(A) 


Ha 


1.3 ± 0.2 


6.6 ± 0.1 


HeH A4686 


0.4 ± 0.1 


6.2 ± 0.1 



Figure 15. As for Figure ImI for Hell A4686. 

In both cases the EW modulation has a maximum at 
phase 0.0 and minimum at phase 0.5. This is 0.25 out of 
phase with respect to the Ha radial velocity curve. In other 
words, the Ha and Hell A4686 EWs peak when the accretion 
disk is between the observer and the secon d ary st ar . 

iGarcia. Callanan. McClintock fc Zhaol il99(il observed 
an Hq EW modulation in GRO J0422-I-32 with a broad min- 
imum at phase 0.5. It peaked when Ha emission line radial 
velocities crossed from the blue to the red side and when 
the secondary star velocities crossed from red to blue. Gar- 
cia et al. attributed the EW changes to an asymmetrically 
emitting disk which may also be asymmetric geometrically. 
The observations were obtained during X-ray outburst and 
decline. 

From comparison to Garcia et al.'s study of GRO 
J0422-I-32 we see the same behaviour in GX 339-4. Hence, 
the EW variations in GX 339-4 may also be due to an asym- 
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metric accretion disk. Investigation of the line profile shape 
over the orbit may help in determining the source of the EW 
variations. This is pursued in Paper II. 

5.3 Full- Width Half Maximum 

The periodograms of the FWHM measured from single- 
Gaussian fits to Ha, Hell A4686, H/3 and the Bowen blend 
showed no significant peaks fsee lBuxtonll20o3) . Hence, the 
FWHM were not considered any further. 

5.4 V/R Ratios 

The periodogram of Ha V/R ratios (Figure I16II shows a 
clear peak at the published orbital period above the 99% 
confidence level. 

The V/R ratios were folded on the orbital period using 
To found in tl5.1l and binned by 0.10 in phase. The results 
are shown in Figure \W\ There is an obvious modulation in 
the ratios. 

Sinusoid curves were fitted to the V/R modulation to 
obtain the semi-amplitude and mean value and their errors. 
The errors on the binned ratios were taken from the standard 
deviation calculated from the least-squares sine fit. 

The minimum Ha V/R ratio is at orbital phase 0.25 
and the maximum at phase 0.75. This is consistent with 
the majority of flux variation originating from the accretion 
disk. The semi-amplitude is 0.11 ± 0.04 and the ratios vary 
around a mean of 1.18 ± 0.03. 

V/R ratios were also measured for Hell A4686 and the 
resulting periodogram is shown in Figure [T7I Although there 
is a single peak at the 99% confidence level it does not cor- 
respond to published orbital period but rather ~ 11.2 hours. 
The V/R ratios of Hell A4686 were folded on this period and 
is presented in Figure [T71 Besides one discrepant data point 
at phase 0.05, the data follow a well-defined sinusoid curve. 
The sold curve is the best-fit sinusoid with a semi-amplitude 
of 0.37 ± 0.09 and mean 0.89 ± 0.06. The maximum falls at 
phase ~ 0.4 and minimum at phase ~ 0.9. 




6 DISCUSSION 

There is no doubt that measuring the orbital period of 
GX 339-4 via disk emission lines will always be difficult. 
This is due to th e low velocity a mplitudes expected from 
a low- inclination dWu et al.ll200J) . This in itself demands 
good quality data, that is, data with high SNR and enough 
spectra to sample the orbital phases well. It was clear in the 
analysis of spectra with low (< 40) SNR and datasets with 
a low number of spectra that one could not obtain satisfac- 
tory periodograms. This is why the results presented here 
have focussed on the dataset of 1998 May 28-31 as it had 
the best data quality. 

In the case of Ha the periodograms of the radial veloc- 
ities, EWs and V/R ratios clearly show an isolated peak at 
the published orbital period at the 99% confidence level. 

For Hell A4686 the radial velocity periodogram did not 
show a single peak but rather a broad double-peak around 
the published orbital period. When one folds the radial ve- 
locities on 14.86 hours anyway one can see a clear devia- 
tion from an otherwise sinusoid curve at phases 0.75-0.95. 



Figure 16. Top: Periodogram of V/R ratios for Ha. y-axis is 
the normalised, reduced x^- Long dashed line shows 99% confi- 
dence level level and dotted line the 68% confidence level. Bottom: 
Folded V/R ratios on the 14.86 hours orbital period. The sine 
curve is best fit to data. See text for fit parameters. 



Whether this is due to some data defect or to a real phys- 
ical source in the binary cannot be discerned at this point 
but it would be interesting to see if it is present in future 
observations. What is clear is that the EW of Hell A4686 
is modulated on the published orbital period and the mini- 
mum and maximum occur at the same phase as the Ha EW. 
Since we expect Ha to originate from the outer disk regions 
and Hell A4686 further in, the source of the EW modulation 
must be the entire disk or at least a large part of it. We leave 
this point here but return to this issue in Paper II when we 
study further the sources of the line profile variations. The 
V/R ratios of Hell A4686 produced a periodogram with a 
single, isolated peak but not at the published orbital period. 
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Figure 17. As for FigureC^for Hell A4686. V/R ratios folded 
on 11.2 hours corresponding to highest peak in periodogram. 



rather at 11.2 hours. It is not obvious why the Hell A4686 
V/R ratios would modulate on a different period than Ha. 
A clue as to why this may be the case is in the asymme- 
try of the line profiles. Ha shows more asymmetry on the 
blue side than Hell A4686. The larger degree of asymmetry 
in Ha would translate into larger amplitudes of V/R ratios 
perhaps making it easier to discern an orbital modulation. 
The V/R modulation of HeH A4686 is fairly clear, however. 

By comparing the phases of minima and maxima of the 
radial velocity, EW and V/R ratio modulations it seems that 
the variations arise from the accretion disk. In Paper II we 
attempt to answer whether these variations are due solely 
to a hotspot, to a disk which is asymmetric geometrically 
and/or in its emitting regions, or both. 



7 SUMMARY AND CONCLUSIONS 

Optical spectra of GX 339-4 was obtained when the system 
was in a high (1998) and low (1999) X-ray state. 

Using interstellar absorption lines from the spectra a 
distance of 4 ± 1 kpc was determined. E[B — V) was mea- 
sured to be 1.1 ± 0.2. These are consistent with previous 
studies. 

Optical colours reddened from 1998 to 1999. This is 
consistent with the accretion disk contributing less blue flux 
to the system due to the decrease in soft X-rays available to 
be reprocessed in the disk. 

As the soft X-rays decreased in the low state the EW of 
Ha, H/3 and Hell A4686 also decreased. Without simultane- 
ous photometry it is impossible to say whether this is mainly 
due to the continuum or line flux variations although the line 
flux must be declining faster than the continuum during the 
low state. 

The FWHM of the Ha, H/3, Hell A4686 and the Bowen 
blend emission lines increased from 1998 to 1999. This sug- 
gests that the line emitting regions moved closer to the com- 
pact object in 1999. The inner regions of the accretion disk 
may have started to cool sufficiently for Balmer emission to 
occur there in the low state. The higher ionisation lines of 
Hell A4686 and the Bowen blend may only obtain sufficient 
photoionising X-rays closer to the compact object. 

The H/3 emission line exhibited redshifted absorption at 
^ A4880 A during both the high and low X-ray states. The 
origin of this absorption feature is not known. 

Emission line profiles vary on the timescales of hours. 
Analysis of the 1998 May high resolution data show that the 
radial velocities, EWs and V/R ratios of Ha modulated on 
the published orbital period. This is the first time that the 
published orbital period has been confirmed by spectroscopic 
data since the study of CCHT92. The EW of Hell A4686 also 
varied on the orbital period but the radial velocities and 
V/R ratios were inconclusive. The radial velocities of the 
Bowen blend modulated on the orbital period but peaked 
0.25 out of phase with respect to Ha. This result may not 
be reliable due to other significant peaks in the periodogram. 

The amplitude of the radial velocity modulation for Ha 
wa s 14 ± 3 km s~^ wh ich is much less than that found 
bv lCallanan et all il992ri . This, in turn, decreases the mass 
function. Taking Porb = 14.86 hours, M2 = 0. 7 Mp, and 
i ^ 7 0° we find that Mi < 40 M©. If i = 15° (IWu et alJ 
I2OOII) then Ml = 5 Mq which places the compact object 
in GX 339-4 into the black-hole category. Determining the 
semi-amplitude velocity from emission lines may not be an 
accurate measur ement of Ki as has been well established 
from CV studies jWaddl 19851) . In order to establish whether 
our measurement of Ki is reasonable future optical spectro- 
scopic studies of GX 339-4 should endeavour to measure the 
semi-amplitude velocities of the emission lines using data of 
high resolution and SNR. The semi-amplitude velocities of 
the same species (e.g. Balmer series) should be measured to 
see if they agree. If not then we need to investigate the phys- 
ical processes which cause such a disagreement which may 
lead to better methods of measuring an accurate Ki from 
emission lines. Until such time, the black-hole categorisation 
of GX 339-4 needs to be taken with caution. 

We may estimate the Roche lobe size (Rl) of the sec- 
ondary star for Port = 14.86 hours. First, let us assume that 
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the secondary is a K dwarf with mass M2 = 0.7 Mq. This is 
a reasonable assumption since the secondary stars in other 
low-mass X-ray binaries with similar periods are of this spec- 
tral type (e.g. A0620-00, Murdin_ct_al. 198C|). If Mi = 5 Mq 
then Rl = 1.3 Kq. If the secondary star in GX 339-4 is a K 
dwarf then it must be slightly evolved to fill its Roche lobe. 
If Ml = 1.4 M0 then Rl = 1.7 Kq. Again, if the secondary 
is a K dwarf, it must be somewhat evolved. 

Future optical spectroscopic observations of GX 339-4 
should endeavor to obtain simultaneous optical photometry 
to ascertain the cause of the EW modulation. The absorp- 
tion trough observed near H/3 should be investigated further 
and it should be resolved whether Hq does or does not ex- 
hibit such a feature. Most importantly, we must obtain as 
many spectra as possible with high SNR and good temporal 
coverage during various X-ray phases and at as many epochs 
as possible. This would also enable us to ascertain whether 
the orbital modulation is observable only if certain physi- 
cal conditions exist in the binary, for example, soft X-ray 
irradiation of a vertically extended hotspot. 

GX 339-4 is a difficult object to decipher due to the lack 
of information about the optical star and the complex be- 
haviour of data. We should, however, continue in our efforts 
to study GX 339-4. This object transitions into all X-ray 
states more often than most and, therefore, provides many 
opportunities to learn about the outburst mechanisms and 
the structure of accretion disks. 
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